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We investigate the superconducting behaviour of Bi doped Pb. Pure lead shows type-I behaviour 
entering an intermediate state in a magnetic field. High dopings of Bi (>3%) lead to type-II be- 
haviour showing a mixed state, where the magnetic field penetrates the superconductor in the form 
of a flux lattice. At intermediate doping, the sample shows both type-I and type-II behaviour de- 
pending on the temperature. This arises because the Ginzburg-Landau parameter k passes through 
its critical value of 1/ \/2 with temperature. 

PACS numbers: 74.25.Dw, 74.55.+h, 76.75.+i, 61.12.Ex 



I. INTRODUCTION 



When a superconductor is brought into a magnetic 
field, its behaviour depends on the balance of its funda- 
mental length scales. The first of these, the penetration 
depth A, describes the length over which the applied field 
penetrates into the superconductor. The second, the co- 
herence length £, roughly gives the spatial extent of the 
Cooper pair. Phcnomcnologically, £ marks the length 
scale over which the macroscopic wave function of the 
supercarriers varies. In the treatment of superconduc- 
tors by Ginzburg and Landau, this wavefunction takes 
the role of the order parameter in the Landau theory of 
phase transitions. The ratio of the two length scales, 
K = A/£, determines the response of the superconductor 
to an applied magnetic field. In the case that k < 1/V2, 
the superconductor is said to be of type-I. Such a type- 
I superconductor will expel a magnetic field completely, 
apart from a layer at the surface of thickness A. This is 
true as long as the field is smaller than some critical field 
B c , associated with the de-pairing current. For applied 
fields higher than B c , superconductivity is lost and the 
field penetrates a sample like a normal metal. The region 
of the thermodynamic B-T phase diagram below B c , in 
which the superconductor behaves as an ideal diamagnet 
is also called the Meissner-state. If the sample geometry 
is such that there are sizable demagnetising effects, even a 
small applied field will exceed the critical field the edges 
of the sample 1 . At these places the field will be able 
to penetrate the superconductor. As a result, normal 
conducting regions with a penetrating field of B c , and 
superconducting regions in the Meissner state (expelling 
the field) will coexist in the sample. This state is called 
the intermediate state. To minimise the surface energy of 
these coexisting domains, the normal state-patches form 
an irregular pattern, similar to that of Weiss-domains 



found in fcrromagnets. 

In addition to the Meissner-state, a type-II supercon- 
ductor (with k > l/v2) shows a second superconducting 
phase, in which it is energetically more favourable for the 
field to penetrate the sample in the form of quantised vor- 
tices of magnetic flux. These flux lines mutually repell 
and form a regular lattice, that is usually of hexagonal 
symmetry. This second superconducting phase consists 
of the mixed state, where the whole of the sample makes 
up a vortex lattice. The mixed state extends up to a field 
B C2 , at which the normalconducting cores of the vortices 
start to overlap. For low fields however, corresponding to 
the Meissner-state, and big demagnetising fields, similar 
arguments to those given above apply. Again, the field 
will exceed the critical field B Cl , above which the sample 
enters the mixed state, at some places. Thus in this situ- 
ation, there will be a coexistence between regions of the 
sample in the Meissner-state, where the field is expelled 
completely, with regions in the mixed state, where the 
field penetrates the sample in the form of a vortex lattice 
with a mean field corresponding to B Cl . 

In the present work, we study samples that have val- 
ues of k close to the boundary between type-I and type- 
II superconductors. By doping Pb samples with small 
amounts of Bi, impurities are introduced to the electronic 
structure. Such impurities lead to a reduction of the 
mean free path of electrons and thus reduce the value of 
£. Therefore small dopings of Bi lead to an increase in 
k. Hence PbBi alloys with sizable amounts (> 3%) of Bi 
are type-II superconductors, whereas pure Pb is a typi- 
cal example of a type-I superconductor. With a suitable 
choice of Bi doping, the coherence length and penetration 
depth are changed in such a way as to permit a crossing 
of the temperature dependences of the critical fields B c 
and B C2 . These critical fields depend solely on the values 
of the penetration depth and the coherence length and 
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are thus strongly affected by doping of the samples. The 
thermodynamic critical field can be calculated to be: 



Be(T) 



$0 



2V2ttA(T)5(T) ' 



(1) 



where $o = h/2e is the flux quantum. Its temperature 
dependence can be calculated in the framework of the 
phenomcnological two- fluid model to be 3 : 



B C (T) = B C (0)(1-(T/T C ) 2 



(2) 



Together with the temperature dependence of A that 
can be independently calculated in the two fluid model 
(A- 2 (T) = A~ 2 (0)(1 - (T/T c ) 4 ), we obtain the tempera- 
ture dependence of £, which determines that of B C2 , given 

by 
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We thus obtain for the ratio of the two critical fields: 
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where k = A(0)/£(0). From this we may then calculate 
a crossover temperature T 7 ™ 77 between type-I and type- 
II behaviour by equating B c and B C2 . This means that 
at temperatures above T 7-77 , the thermodynamic phase 
transition will be at B c and the superconductor will bc 
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FIG. 1. The thermodynamic phase diagram for a thin plate 
of a superconductor with Ko =1. The different transition 
lines corresponding to type-I and type-II superconducting be- 
haviour meet in a multicritical point at a temperature T 1 ^ 11 , 
where the order of the transition changes from first to sec- 
ond order. The phase lines are calculated from the two-fluid 
model (see text). 

of type-I. Below T 7 " 77 however, the phase transition to 
the normal state happens at B C2 and thus the supercon- 
ductor is of type-II. From the above arguments, we obtain 
for T 1 - 
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T 1 - 11 = T c (V2k - l) 1 
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which only has a solution below T c for 1 / V2 < Ko < \/2, 
putting a severe limit on the observability of the effect. 
More exact values may be determined from microscopic 
theories such as those of Hclfand and Werthamer, taking 
the electron mean- free paths into account 4 . The results 
however do not differ greatly and would only be available 
numerically, thus obscuring the nature of the effects. In 
the B-T phase diagram, this point at which the funda- 
mental behaviour of the superconductor changes may be 
interpreted as a multicritical point. An overview of this 
situation can be seen in Fig. 1, where we show the B- 
T phase diagram of a superconductor with Ko =1. The 
different transition lines B c , B C2 and B Cl cross at the 
multicritical point at T 7 ~ 77 . Furthermore, the phase 
transition at B c in a type-I superconductor is of first 
order, whereas that into a type-II superconductor at B C2 
is of second order 3 . The third phase line B Cl , separat- 
ing the mixed state from the intermediate-mixed state in 
type-II superconductors depends on the value of k and 
the nature of the interactions between flux lines. For 
low k and attractive vortex interactions, the transition 
is thought to be of first order 5 , otherwise it is second 
order. Thus we expect to observe a small region of coex- 
istence between type-I and type-II superconducting be- 
haviour around T 7 ~ 77 . We will discuss this transition 
from type-I to type-II behaviour in more detail below in 
the discussion of the experimental findings. 



II. EXPERIMENT 



We investigated four different samples of PbBi alloys, 
with different doping levels of Bi. For a typical example 
of type-I behaviour, we used a pure Pb sample and for a 
determination of the type-II properties we investigated a 
sample with a doping of 5% Bi. The samples with inter- 
mediate doping, 1.25 and 1.5%, showed a transition from 
type-I to type-II behaviour at different fields and tem- 
peratures. The samples all were polycrystalline platelets 
of approximate dimensions 40x25x1 mm 3 , and an ellip- 
soidal cross-section. The samples were characterised with 
a vibrating sample magnetometer (VSM) and showed a 
superconducting transition, T c of 7.2 K in zero applied 
field. 

The muon spin rotation (^(SR) experiments were car- 
ried out at the MUSR instrument of the ISIS facility at 
the Rutherford Applcton Laboratory (RAL), UK. The 
ISIS facility provides a pulsed beam of positive muons 
at a frequency of 50 Hz with a pulse-width of ~ 70 ns. 
These muons are created in the decay of pions at rest and 
hence are almost fully spin polarised, antiparallel to their 
momentum, and have a kinetic energy of ~ 4 MeV. The 
neutron scattering measurements were done at the Insti- 
tut Laue-Langevin (ILL), France, using the instrument 
Dll. ILL provides the world's most intense source of 
cold neutrons, where in our investigations we have used 
neutrons of a wavelength of 1.9 nm. For the very large 
structures in the intermediate mixed state, we needed to 
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FIG. 2. a) The local field distribution in the intermediate 
state. A background of muons stopping in the cryostat shows 
an applied field of 20 mT, while the muons in the sample pre- 
cess at a frequency corresponding to 38 mT. This corresponds 
to the critical field B c at the temperature of 5.3 K, where the 
measurement was carried out. b) The local field distribution 
in the intermediate mixed state. The applied field can be 
seen to be 40 mT. In addition, there is a wide distribution of 
local fields corresponding to a vortex lattice at higher fields. 
The mean field of this flux lattice corresponds to the critical 
field at the transition from type-I to type-II behaviour T 1 ^ 11 , 
if the overall field expulsion is taken into account (see text). 
The measurement was carried out at a temperature of 4.3 
K. c) The field distribution for a vortex lattice, in the mixed 
state. The field distribution has a very characteristic shape, 
with a tail extending to fields higher than the average. The 
measurement was taken in an applied field of 30 mT in the 
type-II sample containing 5 % Bi. 

observe small-angle scattering, for which Dll with a dis- 
tance between the sample and the detector of up to 40 m 
is ideal. In the investigations described below, we have 
used a collimation length of 20.5 m and a sample-detector 
distance of 20 m. No beam stop was used. The neutrons 
are detected with a position-sensitive detector of 0.8x0.8 
m 2 area and pixel size 10 mm. In both experiments, the 
field was applied at 45° to the shortest direction of the 



sample, in order to have sizable demagnetisation fields. 

In a transverse- field fiSK experiment, positive muons 
come to rest in the bulk of a sample, where the field is 
applied perpendicular to the initial spin polarisation of 
the muons. Due to their magnetic moment, the muons 
then perform a Larmor-precession with a frequency pro- 
portional to the applied field. After a mean life-time of 
2.2 f/,s, the muons decay into a positron (and two neu- 
trinos). Due to parity violation in the weak decay of 
the muon, the decay-positron is preferentially emitted in 
the direction of the muon spin. Therefore the number of 
positrons observed at a fixed position, in the precession 
plane of the muons, will oscillate with time reflecting the 
Larmor-precession of the muons, in addition to an expo- 
nential decay, due to the radioactive decay of the muons. 
For a distribution of fields over the sample, as is the case 
for a vortex lattice, the oscillations due to the different 
fields superimpose giving rise to a depolarisation of the 
muons. In the case of large field gradients, the depolari- 
sation rate may be fast enough to reduce the initial polar- 
isation outside of the experimental time window. In the 
mixed state of a superconductor with short penetration 
depth, this may strongly reduce the signal observed by 
/iSR. By studying the Fourier-transform of the number of 
decay positrons with time, we may however still observe 
the probability distribution of the internal fields, assum- 
ing a random distribution of muons over the sample. In 
order to avoid spurious noise in the Fourier-transform 
arising from the finite time window of observation, we 
use a maximum-entropy algorithm in our analysis 9 . This 
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FIG. 3. The thermodynamical critical field B c , as deter- 
mined from /iSR in the intermediate state. The measure- 
ments are for all three samples exhibiting type-I behaviour. 
All curves were obtained with several applied fields, where 
the different determinations overlap (see text). It can also be 
seen that B c is the same for all samples in spite of the change 
in ft. 

also prevents statistical noise at long times, due to low 
count rates, from spreading over the whole frequency do- 
main. This is mainly because our algorithm never ac- 
tually Fourier-transforms the raw data, but rather solves 
the inverse problem of calculating the frequency distribu- 
tion and compares it to the time-data. A more detailed 
discussion can be found in Ref. 8 . 



For the vortex lattice, in the mixed state, the field 
distribution has a very specific shape, with a tail extend- 
ing to fields higher than the average 10 . This is due to 
the high fields in the vortex cores. In the intermediate 
state, we expect a single-frequency signal at the critical 
field B C (T). The signal size in the intermediate state is 
given by the area fraction of the normal state domains. 
Neglecting the effects flux expulsion this is simply pro- 
portional to B / B C (T), leading to a decrease in signal 
with decreasing temperature. Demagnetising effects due 
to the experimental geometry of having the field at 45° 
to the plate however lead to a much stronger reduction of 
signal size with temperature. This leads to a dependence 
on the critical field cx 1/(B—B C ), in reasonable agreement 
with the experimental findings. This may also be ob- 
served from the asymmetry in decay positrons upstream 
and downstream of the beamline. This corresponds to 
muons not precessing in a field that hence stopped in a 
Meissner domain in the sample. The temperature depen- 
dence of this signal is complementary to that observed 
from the normal state domains. 

In a neutron scattering experiment, the neutrons are 
incident almost parallel to the applied field. The align- 
ment of the beam with respect to the field was done by 
observing the diffraction pattern of the vortex lattice in 
a sample of Nb to an accuracy of 0.1° both vertically and 
horizontally. Due to its magnetic moment, the neutron 
is scattered by gradients of magnetic field 2 . Therefore 
we may observe the domain pattern in the intermediate 
state in terms of a correlation function. In general, the 
scattered neutron intensity is given by 7 



given below. This structure factor can be calculated to 
be 11 
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where Q is the scattering vector, R(Q) is a resolution 
function, describing the distribution of neutrons around 
the nominal scattering vector Q. Finally, da/dQ is given 
by 7 : 
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dn 
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where F{Q) is a form factor describing the scattering 
of a basic building block of the structure and S(Q) is 
a structure factor, describing the large scale structure. 
The fourier transform of the structure factor presents a 
measure of the 'pair correlation' or distance distribution 
function in the large scale structure. Thus eq. 7 corre- 
sponds to a convolution of the correlation functions of 
the two parts of the scattering. 

In the intermediate state, this is a similar problem to 
that encountered in the study of vesicles of amphiphilic 
molecules that may build large scale structures. In con- 
trast to these amphiphilic molecules however, the micro- 
scopic structure of the magnetic field domains is very 
simple. This microscopic structure is described by the 
London equations describing the magnetic behaviour of 
superconductors. Thus in the general expression for the 
scattered intensity, the form factor will be given by the 
well known London form factor. For the structure factor, 
we have used a model of randomly oriented chains with a 
radius of gyration R g in the analysis of the measurements 



S C (Q) = 



2(exp(-R 2 g Q 2 ) - 1 + R 2 g Q 2 



(8) 



In the limit of large R g Q 1 this form factor has the asymp- 
totic dependence S C (Q) cx Q~ 2 , corresponding to a ran- 
dom walk arrangement of the chains. From the determi- 
nation of R g , we may then determine the contour length 
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FIG. 4. a) The scattered neutron intensity as a function of 
the scattering vector Q. A diffraction pattern above T c has 
been subtracted. The solid line presents a fit to a model using 
the structure factor of eq. 8. Thus we determine the radius of 
gyration and the interface area of the domain structure (see 
text), b) These measurements were done in the intermedi- 
ate-mixed state. In addition to the structure factor from the 
domains in the intermediate-state, there is a strong peak cor- 
responding to Bragg diffraction from the vortex lattice also 
present. The lattice spacing corresponds to a field of 45 mT, 
the critical field at T /_/I ,the transition from type-I to type-II 
behaviour. The solid line represents a fit to the same model 
as above (see text) . 

of a chain, from L = 6R 2 /b, where we assume the statis- 
tical segment length b to be the penetration depth. This 
model will break down below the length scale, where the 
magnetic domains become stiff. On this scale, a model 
of randomly distributed, solid rods is more appropriate. 



A 



In this case, the structure factor has the asymptotic de- 
pendence S r oc Q^ 1 . For the magnetic domains, in the 
intermediate state, this length scale will be of the order 
of the penetration depth (A ~ 90nm) . This makes it nec- 
essary to observe very small scattering vectors, in order 
to resolve the domain structure in the intermediate state. 
Due to the very small scattering angles the observed scat- 
tering distribution is also influenced by the width of the 
main beam. This may be taken into account by fitting 
a Gaussian distribution to the main beam in the normal 
state. The intensity of the main beam in the supercon- 
ducting state is then reduced by those neutrons scattered 
according to S C (Q), whereas we assume the width of the 
main beam to be that in the normal state. From the 
radius of gyration thus determined we may then obtain 
the contour length of a domain in the intermediate state 
from 



R 2_Lb 
9 ~ ~6~' 



(9) 



where the Kuhn segment b is of the order of the pen- 
etration depth A and L is the contour length. In Lan- 
dau's treatment of the intermediate state 12 , the separa- 
tion of the different domains in an equilibrium situation 
depends only on the thickness of the plate, the surface 
energy length 6, usually associated with the coherence 
length in type-I superconductors and a universal function 
f(b), where b is the reduced field B/B c . This function is 
calculated from the different energy scales of the order 
parameter and the magnetic field and is given by 12 
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-(1 + b 2 yin(l + b A ) - 4b 2 ln(8b)}. 

The separation of the different domains in the interme- 
diate state is then given by the expression 
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where d is the thickness of the plate. In the case of in- 
terest here, where k is close to the boundary of 1/V2, 
the surface energy length becomes very small, such that 
the separation of the domains may of similar order as 
the separation of vortices in the mixed state and thus 
observable in neutron scattering experiments. This will 
be discussed further in the context of the experimental 
findings. 

In the mixed state, the situation is simpler, as the 
scattering now is from a lattice. Therefore the struc- 
ture factor will exhibit narrow peaks when the scattering 
vectors equal a reciprocal lattice vector. The big lat- 
tice spacing in the vortex lattice (<~ 200 nm at 60 mT) 
leads to very small angles of scattering (~ 1°). In the 
intermediate-mixed state, there will be both structures, 
that of a vortex lattice and the characteristic domain 
structure of the intermediate state. The structure factor 
will thus be made up of a combination of both Bragg and 
domain-scattering. 



III. RESULTS AND DISCUSSION 



A typical /iSR measurement of the intermediate state 
in a type-I superconductor is shown in Fig. 2a) for the 
case of the 1.25 % Bi sample at a temperature of 5.3 
K. Due to a small background of muons stopping in the 
cryostat, we observe the applied magnetic field (20 mT). 
There is also a second peak at the critical field B c corre- 
sponding to the temperature of 5.3 K. This field can be 
seen to bc 38 mT in the figure. Thus we can map the 
phase boundary B C (T) in the type-I state, by measuring 
the local field in the sample at different temperatures. 
This is shown in Fig. 3 for the pure Pb sample, as well 
as the 1.25% Bi and 1.5% Bi samples in the type-I state. 
Due to the loss of signal with increasing field gradient, the 
measurements were performed over an overlapping range 
of applied fields. The observed B C (T) curve is found to be 
independent of the applied field. Furthermore, the tem- 
perature dependence of B c is in good agreement with the 
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FIG. 5. The distribution of local fields as observed by /iSR 
in the 1.25% Bi sample at a temperature of 4.7 K. The applied 
field is 40 mT. Moreover, one observes a sharp peak at a field 
corresponding to the critical field at that temperature (see 
Fig. 3) . Finally one observes a very wide distribution of fields 
below B c . This corresponds to flux lines that have nucleated 
out of the normalconducting regions in the intermediate state. 
Therefore in this situation, one observes a coexistence of the 
intermediate and the intermediate-mixed state, as would be 
expected at a multicritical point (see text). 

expectations from the two-fluid model, as can be seen 
from the solid line in Fig. 3, representing a fit to eq. 2. 
Moreover, it can be seen in the figure that for the doped 
samples the value of B c remains unaltered. Consider- 
ing the dependence of B c on the values of the coherence 
length and the penetration depth (eq. 1), we note that 
X'C = A£, where A' and £' are the values of the pene- 
tration depth and coherence length after doping. This 
agrees with the expectations of microscopic theories on 
the influence of impurities on the superconducting pa- 
rameters in the temperature region of interest 4 . 

A typical SANS measurement of the intermediate state 
can be seen in Fig. 4a). This figure shows the 1.25% Bi 
sample at 5.3 K, as Fig. 2a) does for the ^SR measure- 
ments. The intensity shown presents a subtraction of 



the scattering observed at 5.3 K and that above T c . The 
data presents a radial average over the detector, giving 
the dependence of the scattered intensity on the abso- 
lute value of Q. The solid line in Fig. 4, represents a fit 
to eq. 7, where the structure factor was assumed to be 
given by eq. 8. Furthermore we have taken into account 
the effects arising from the very small scattering angles 
as discussed above. The value for the penetration depth 
used in calculating the London form factor was obtained 
from the /iSR measurements on the vortex lattice in this 
sample. From this fit we determine the scattered inten- 
sity and the radius of gyration of the domain structure. 
Here, the scattered intensity presents a measure of the 
interface area between the domains and the radius of gy- 
ration roughly gives the size of the domains. This shows 
the very complementary results obtained from /iSR and 
SANS. 

Fig. 2b) shows the field probability distribution ob- 
served in /nSR in the intermediate-mixed state. The mea- 
surement was done in the 1.25% Bi sample at a tempera- 
ture of 4.3 K with an applied field of 40 mT. In addition to 
the background of muons stopping in the cryostat, we ob- 
serve a signal arising from a flux line lattice with a mean 
field much higher than the applied. We will see below 
that the internal field corresponds to £? C (T 7 ~ /7 ), the crit- 
ical field at the temperature where the superconducting 
behaviour changes. The same behaviour can be seen in 
Fig. 4b), where the scattered neutron intensity is shown 
at the same temperature, in an applied field of 35 mT. 
As in Fig. 4a), the scattered intensity above T c has been 
subtracted and the London form factor has been taken 
into account. Similar to Fig. 4a), we observe scattering 
from the domain structure are low Q, that is similar to 
that observed in Fig. 4a), but in addition there is a peak 
in scattered intensity arising from the vortex lattice in 
the intermediate-mixed state. The plane spacing of this 
lattice, as obtained from the Q-value of the peak, corre- 
sponds to a mean field of a triangular lattice of 42 mT, 
in good agreement with the results from /iSR. The plane 
spacing of the vortex lattice is connected to the field by 
the fact that each vortex line carries a quantum of flux. 
Thus the plane-spacing is given by d = (v / 3/2$o/-B) 1 ^ 2 
for a triangular lattice. From a comparison of Fig. 2b) 
with Fig. 2c), where the field distribution of a vortex 
lattice is shown, as obtained in the type-II sample con- 
taining 5 % Bi, it can be seen that the field distribu- 
tion in the intermediate-mixed state does not perfectly 
match that of the mixed state. This may be understood 
from a comparison of the lattice spacing with the sepa- 
ration of the domains in the intermediate state. As will 
be shown below, they are of the same order of magnitude 
(L ~ 10A), which may lead to imperfections in the vortex 
lattice resulting in a more symmetric lineshape. 

The value of the mean internal field in the 
intermediate-mixed state can be understood by consider- 
ing the transition between type-I and type-II supercon- 
ducting behaviour. The penetration of the field in the 
form of a vortex line arises from a competition between 
the surface energy gained by expelling the magnetic field 
to the cost of an inclusion of normalconducting material. 
The first will be roughly proportional to A, whereas the 



second will be proportional to £ (see Ref. 6 ). From this 
argument, it is seen that for a type-I superconductor, it 
is energetically unfavourable to build up domain bound- 
aries. A rigorous calculation of these surface energies 
leads to the criterion of n — 1/ V2 for the crossover from 
type-I to type-II behaviour when these energies are equal. 
We now consider the case of the intermediate- mixed state 
found in our samples at below T 1 ^ 11 . On cooling the 
sample, it is in the intermediate state with normalcon- 
ducting domains of mean field corresponding to B C (T). 
At the temperature T /_// , the ratio of the surface en- 
ergies passes through one. Thus as discussed above the 
formation of normalconducting-superconducting bound- 
aries becomes energetically favourable. This corresponds 
to a sign-change of an effective surface tension of the 
magnetic field in the domains, leading to the nucleation 
of vortices from the normalconducting domains, due to 
flux quantisation in the superconductor. The domain 
structure however, remains the same, which is why the 
mean field of the vortex lattice corresponds to the value 
of B C (T J_JT ), due to flux conservation. 
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FIG. 6. The scattered neutron intensity from the domain 
structure in the intermediate and intermediate-mixed state as 
a function of temperature at three different applied fields in 
the 1.25% Bi sample. The scattered intensity was obtained 
by a fit using the model given by eq. 8. This presents a mea- 
sure of the interface area of the domains, thus describing the 
ruggedness of the flux distribution in the sample (see text). As 
would be expected the scattered intensity shows a maximum 
around the temperature T 1 ^ 11 , where vortices are nucleating 
from the normal domains (see text). 

As the order of the superconducting transition changes 
between a type-I and type-II superconductor in a mag- 
netic field, the point in the B-T phase diagram at T 1 ^ 11 
may also be interpreted as a multicritical point. In this 
scenario, we would expect to see a co-existence of the in- 
termediate and the intermediate- mixed state over a short 
temperature interval. That this is in fact the case can be 
seen in Fig. 5, where we show the local field distribu- 
tion observed from ^SR at a temperature of 4.7 K in the 
1.25% Bi sample. As can be seen in the figure, in addi- 
tion to the sharp peaks at the applied field and B c , there 
is a field distribution at fields just below B c . This coex- 
istence of the intermediate and the intermediate-mixed 
state is observed over a temperature interval of 0.4 K. 



Within our temperature resolution of 0.1 K there is no 
hysteresis observable in the ratio between the two states. 

From the onset of superconducting behaviour at higher 
fields, we determine the value of B C2 , which together with 
the value of B c determined from the Pb sample results 
in a determination of kq- For the two samples showing a 
transition from type-I to type-II behaviour, this results 
in 4' 25 = 1.02(2) and kJ 5 = 1.21(3). Using eq. 5, we 
may then estimate the transition temperature T /_7/ , re- 
sulting in = 4.8(2) K and T^ 11 = 6.1(1) K. These 
values are in excellent agreement with those obtained di- 
rectly in both fiSR and SANS from the appearance of 
the intermediate-mixed state. These directly determined 
values are T^ 1 = 4.7(1) K and T^ 11 = 6.2(1) K. This 
also indicates that the phenomenological two-fluid model 
gives a fair description of the superconducting behaviour 
of these PbBi alloys. 

We now briefly return to the observation of the in- 
termediate state with SANS. As already noted above, 
from the structure factor we determined in the interme- 
diate state, we may obtain a measure of the interface area 
from the total scattered intensity, as well as the radius 
of gyration of the domain structure. Both these quan- 
tities depend upon temperature. At high temperatures, 
close to T c , the critical field is close to the applied field, 
such that most of the sample area will be in the nor- 
mal state and the domain interface will be small. At 
low temperatures, the critical field is much higher than 
the applied field and only a very small area will be in 
the normal state. Thus again the domain interface will 
be small. Therefore the scattered intensity will exhibit a 
maximum with temperature, where the domain structure 
is most convolved. This can be seen in Fig 6, where we 
show the temperature dependence of the scattered inten- 
sity in the 1.25% Bi sample at different applied fields. As 
can be seen in the figure, the intensity shows a marked 
peak at the transition temperature T 1 ^ 11 and then de- 
creases slowly at low temperatures. Such a sharp peak in 
the region of coexistence between the intermediate and 
the intermediate-mixed state is possibly due to the strong 
increase in surface area of the magnetic flux structures 
when the two states coexist. At this point, the domain 
structure does not change, but the normalconducting re- 
gions are nucleating vortex lines, leading to both an en- 
hancement in the surface area and the size of the mag- 
netic domain structures. The width of the coexistence is 
similar to that found by /iSR. The decrease in intensity 
is faster at higher applied fields. This may arise from the 
fact that these high fields are already sufficiently close 
to B Cl , such that the superconductor would simply be in 
the mixed state at low temperature. Hence there would 
be no scattering from domains in that situation. 

In the same way, we may investigate the size of the 
domains from the radius of gyration observed in the in- 
termediate state. Using the relationships between the 
radius of gyration determined from the fits to the scat- 
tered neutron intensities with the separation of the do- 
mains and Landau's expression for this separation, we 
may also determine the universal function f(b). With 
the known thickness of the sample of 1 mm, there are 
two adjustable parameters left, whose order of magni- 



tude is known. These are the Kuhn length of the seg- 
ments of the chains that is roughly given by the pene- 
tration depth and the surface energy length that has to 
be somewhat smaller than the coherence length. Using a 
segment length equal to the penetration depth, the sur- 
face energy length S = 10 nm has to be chosen to obtain 
reasonable agreement with the expectations from Landau 
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FIG. 7. The function f(h) from Landau theory, as deter- 
mined from the radius of gyration of the normalconduct- 
ing domains in the intermediate and the intermediate-mixed 
state. The temperature dependences at different fields allow 
us to determine f(h) over a big range of reduced fields. In or- 
der to obtain reasonable agreement of the size of the domains 
with the predictions of Landau theory, a surface energy length 
of 10 nm, much smaller than the coherence length had to be 
assumed. This is however not unreasonable, given the close- 
ness of k to \j\p2 (see text). 

theory for f(b). These results can be seen in Fig. 7 for 
the 1.25 % sample. The use of the temperature depen- 
dences at the different fields results in a determination 
of f(b) over a big range of values. The value of 5 used 
to obtain the agreement with Landau theory is an order 
of magnitude smaller than the coherence length. This is 
however not unexpected, as the value of k is very close to 
its boundary value of 1/ \pl and hence the surface energy 
length is expected to be small as discussed above. 



IV. CONCLUSIONS 

We have measured the behaviour of type-I and type-II 
superconductors in a magnetic field. It is shown that 
there is a transition from type-I to type-II behaviour 
with temperature. This transition happens at a tem- 
perature well defined in the framework of the two-fluid 
model and can be shown to only arise for superconduc- 
tors with a Ginzburg-Landau parameter in the interval 
1/V2 < k < V2 at low temperature. This transition 
may be seen as a multicritical point in the B-T phase 
diagram of superconductors, as the order of the super- 
conducting transition is changed from first to second or- 
der. Therefore we also observe a coexistence of the dif- 
ferent behaviours around the transition with both SANS 
and /iSR. We furthermore use the complementarities of 
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SANS and ^iSR to study the intermediate state in type-I 
superconductors. Using SANS, we observe the domain 
structure via the interface area given by the scattered in- 
tensity and the domain sizes via the radius of gyration, 
where there is qualitative agreement with the treatment 
of the intermediate state by Landau. The local fields are 
investigated with /iSR, where we can map the critical 
field as a function of temperature. 
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